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ABSTRACT 

O 

^ I We present constraints on several cosmological parameters from a combined 

O I analysis of the most recent Cosmic Microwave Background anisotropy data and 

^ ' the Sloan Digital Sky Survey cluster mass function. We find that the combination 
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of the two data sets breaks several degeneracies among the parameters and pro- 
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vides the following constraints: cxg = 0.76 ± 0.09, Qrn = 0.26+|]:E]^, h = 0.66+°'°^ 



n = 0.96 ± 0.05, = 0.07^°;°^. 
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1. Introduction 



Q ■ The last few years have seen a spectacular increase in the amount and quality of available 

^ cosmological data. The new results on the Cosmic Microwave Background angular power 

ci • spectrum (Netterfield et al. 2001; Halverson et al. 2001; Lee et al. 2001; Pearson et al. 

2002; Scott, et al. 2002; Benoit et al. 2002) have confirmed the theoretical prediction of 
^ ■ acoustic oscillations in the primeval plasma and constrained theories of large-scale structure 
■ formation (see e.g. Wang, Tegmark & Zaldarriaga 2001). At the same time, early data 
from the 2dFGRS (Percival et al. 2001) and the Sloan Digital Sky Survey (York et al. 2000; 
Stoughton et al. 2002) galaxy redshift surveys have provided an unprecedented view of the 
large-scale structure of the universe as traced by galaxies. 

Combined analysis of these independent CMB and galaxy data sets are placing strong 
constraints on some of the fundamental cosmological parameters (Bahcall et al. 1999; Efs- 
tathiou et al. 2001; Lahav et al. 2002; Melchiorri & Silk 2002). Together with the high-redshift 
supernovae results (Perlmutter et al. 1999; Filippenko & Riess 2000), a concordance model 
of a flat, low-density Cold Dark Matter cosmology has become the current paradigm. 
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The goal of these analyses is to determine the precise values of the cosmological param- 
eters of the A-CDM model. Due to 'cosmic degeneracy', the CMB data alone are unable to 
place tight constraints on several fundamental parameters, such as the r.m.s. amplitude of 
the mass fluctuations cig, the Hubble parameter h and the optical depth of the universe Tc, 
even if one restricts the analysis to a flat universe. 

In the present Letter we combine the most recent CMB anisotropics data with the 
constraints obtained from the mass function of clusters of galaxies determined from early 
commissioning imaging data of the SDSS (Bahcall et al. 2003) to break the degeneracy 
among the cosmological models, allowing a determination of the best-fit values for individual 
parameters. 

2. Method 

2.1. CMB data and analysis 

We consider a template of flat, adiabatic, A-CDM CMB spectra computed with CMB- 
FAST (Seljak & Zaldarriaga 1996), samphng the various parameters as follows: the physical 
density in cold dark matter O^cdmh^ = '^cdm — 0.05, ...0.40, in steps of 0.02; the physical den- 
sity in baryons Qft/i^ = ujb — 0.009, 0.030, in steps of 0.003; and the cosmological constant 
Qa = 0.5, 0.95, in steps of 0.05. For each set of these parameters, the scalar spectral index 
n is varied over the relevant inflationary values of n = 0.8, 1.2, in steps of 0.02. The value 
of the Hubble constant is not an independent parameter, since h = {oOcdm + ^b) /(I — ^a); 
we include a top-hat prior h = 0.7 ± 0.2 (Freedman et al. 2001). Only models with age 
to > 11 Gyrs are considered. 

We allow for a reionization of the intergalactic medium by varying the Compton optical 
depth parameter t^. in the range Tc = 0.0, ...,0.45 in steps of 0.05. High values of Tc are in 
disagreement with recent estimates of the redshift of reionization 2:^6 ~ 6 ± 1 (see e.g. Fan et 
al. 2001; Gnedin 2001), which point towards Tc ~ 0.05 — 0.10. However, since the mechanism 
of reionization is still not clear, we allow this parameter to vary freely within the above 
conservative range. As shown below, the combination of the CMB and CMF data provides 
an independent constraint on this parameter. 

For the CMB data, we use the recent results from the BOOMERanG-98 (Netterfleld et 
al. 2001), DASI (Halverson et al. 2001), MAXIMA-1 (Lee et al. 2001), CBI (Pearson et al. 
2002), VSA (Scott, et al. 2002) and Archeops (Benoit et al. 2002) experiments. 

The power spectra from these experiments were estimated in 19, 9, 13, 14, 10 and 16 
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bins respectively. For the CBI, we use the data from the MOSAIC configuration (Pearson 
et al. 2002), spanning the range 2 < ^ < 1500. We also use the COBE data from the 
RADPACK compilation (Bond, Jaffe & Knox 2000). 

For the Archeops, CBI, DASI, MAXIMA-I and VSA experiments we use the pubhcly 
available correlation matrices and window functions. For the BOOMERanC experiment 
we assign a flat interpolation for the spectrum in each bin £(£ + l)C^/27r = Cb, and we 
approximate the signal Cb inside the bin to be a Gaussian variable. The CMB likelihood 
for a given theoretical model is defined by —2\\iC^^^ = {C%' — C^)Mbb'{C*^i — C^,) where 
Mbb' is the Gaussian curvature of the likelihood matrix at the peak. 

We consider 7%, 10%, 4%, 5%, 3.5% and 5% Gaussian distributed calibration errors (in 
AT) for the Archeops, BOOMERanG-98, DASI, MAXlMA-1, VSA, and CBI experiments 
respectively and we include the beam uncertainties by the analytical marginalization method 
presented in Bridle, et al. (2001). Finally, we rescale the spectrum by a prefactor Cio, 
assumed to be a free parameter, in units of C^^^^ . 

In order to constrain a parameter x we marginalize over the values of the other param- 
eters y. This yields the marginalized likelihood distribution 



The central values and Icr limits are then found from the 16%, 50% and 84% integrals of 



We use the cluster mass function (CMF) obtained from the early SDSS commissioning 
data (Bahcall et al. 2003). This CMF was derived from 294 clusters in the redshift range z 
= 0.1-0.2 selected by the Hybrid Matched Filter (HMF) method. The HMF mass function 
was compared with large scale cosmological simulations as well as with model predictions, 
as discussed in Bahcall et al. (2003). If the Hubble constant and spectral index are kept 
constant at h=0.72 and n=l, then the best-fit relation between amplitude and density can 
be summarized as crgQ'^ = 0.33. The 68, 95, and 99% confidence contours (allowing h and 
n to vary as in the previous section) are shown by the dashed curves in Figure 1. 

The shape of the cluster mass function only partially breaks the degeneracy between flrn 
and (Tg in the above relation: low values of as yield a steeper CMF shape at the high-mass 
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2.2. Cluster Mass Function Analysis 
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end (i.e., fewer high mass clusters) than do low Qm values (which produce a flatter CMF 
shape). The cluster mass function prefers a low value for the mass density parameter and 
a relatively high value for the amplitude as'- the best-fit parameters for the HMF clusters 
are ^2^=0. 18 and o"8=0.92; similar results are obtained for SDSS clusters selected by the 
the maxBCG method (see Bahcall et al. 2003). The above ag — Qm relation is consistent 
with recent results from the X-ray cluster temperature function (see Seljak 2001), and from 
cosmic shear lensing observations (Ikebe et al. 2002; Reiprich & Bohringer 2002; Hamana et 
al. 2002; Jarvis et al. 2002). It also implies that for a mass density of Q^=0.3 the relevant 
amplitude is erg = 0.7. 

The observed CMB spectrum of fluctuations suggests a lower amplitude value for erg 
(which is, however, degenerate with the optical depth parameter for CMB), and a somewhat 
larger value for the mass density parameter (Q^n ~ 0.3 for h—0.72), than given by the clus- 
ter mass function above (see e.g. Lahav et al. 2002; Melchiorri & Silk 2002, and references 
therein). However, the CMB and the cluster mass function results are consistent with each 
other within one sigma. Combining the CMB and CMF data will clearly result in interme- 
diate values for the cosmological parameters, shifting the CMF constraints presented above 
towards a somewhat lower amplitude and higher mass density regime. 

We combine the cluster mass function results with those of the CMB by multiplying 
the two likelihoods C^mb^cmf^ 

using the same range of parameters {ftm, <^s, h, n) and 
marginalizing over the nuisance parameters as discussed in the previous section. 



3. Results 

The main result of our analysis is presented in Figure 1, where we plot likelihood contours 

in the Vtm — erg plane for the two datasets, separately and combined. It can be seen that 
both the CMB and CMF datasets are affected by degeneracies between and Q„i- In the 
case of the CMF, an increase of Q„i results in a larger number of clusters, and must be 
reduced to bring the predicted CMF back in line with observations. On the other hand, the 
CMB is only weakly sensitive to the tradeoff between density and amplitude, so models with 
higher and as can be in agreement with CMB data. The degeneracy in the CMB dataset 
is therefore opposite to the one in the CMF data, and the two measurements complement 
each other. The combination of the two datasets provides the constraints: — 0.261;q;o7 
and (78 = 0.76 ± 0.09 at la confldence level, as shown in Figure 1. 

The CMB-I-CMF combination can also break additional degeneracies. The optical depth 
Tc and the Hubble parameter h are better determined after the inclusion of the CMF data. 
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Prom the combined analysis, we obtain (at 1 — cr C.L.) h = 0.66to'^l, n = 0.96 ± 0.05 and 
Tc = 0.07+[J:[!^; these values can be compared with h = O.rOtoH, n = 0.98l^;g^ and < 0.29 
from the CMB-only analysis. 

4. Conclusions 

We combine the Cosmic Microwave Background anisotropy and SDSS cluster mass func- 
tion data to produce constraints on several cosmological parameters. The complementary 
nature of the two data sets breaks existing degeneracies among cosmological parameters. 
The CMB data tend to indicate a higher and a lower amplitude cTg than suggested by 
the cluster data; thus the combined result for is pulled to a lower value than suggested by 
the CMB alone, and the amplitude cg is lower than suggested by the cluster data alone. We 
find the combined data suggest a mass density = 0.261o;q7; a normalization of the matter 
power spectrum erg = 0.76 ± 0.09; Hubble parameter h — 0.66to;o6; a nearly scale-invariant 
spectrum of primordial fluctuations n — 0.96 ± 0.05; and an optical depth of the universe 
Tc = 0.071q;q5. We have restricted the analysis to flat universes; if this is the case, the density 
of mass-energy in the universe is dominated by a form other than ordinary matter. 

We acknowledge a grant from the Oxford- Princeton Research Partnership Initiative. AM 
is supported by PPARC; PB is supported by the National Computational Science AUiance 
under NSF Cooperative Agreement ASC97-40300. 
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Fig. 1. — Confidence contours in the Qm — <^8 plane from tlie CMB+CMF combined analysis. 
Bold dashed and faint dashed curves represent the 68, 95, and 99% confidence contours from 
the CMF and CMB data, respectively. The shaded regions represent the corresponding 
ranges allowed by the combined analysis. 



